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ABSTRACT 

A powerful way to observe directly the solid state inventory of dense molecular clouds is by infrared 
spectroscopy of background stars. We present Spitzer/IRS 5-20 /im spectra of ices toward stars behind 
the Serpens and Taurus molecular clouds, probing visual extinctions of 10-34 mag. These data provide 
the first complete inventory of solid-state material in dense clouds before star formation begins. The 
spectra show prominent 6.0 and 6.85 /im bands. In contrast to some young stellar objects (YSOs), 
most (~75%) of the 6.0 /im band is explained by the bending mode of pure H2O ice. In realistic 
mixtures this number increases to 85%, because the peak strength of the H2O bending mode is very 
sensitive to the molecular environment. The strength of the 6.85 /im band is comparable to what 
is observed toward YSOs. Thus, the production of the carrier of this band does not depend on the 
energetic input of a nearby source. The spectra show large abundances of CO and CO2 (20-40% with 
respect to H2O ice). Compared to YSOs, the band profile of the 15 /im CO2 bending mode lacks 
the signatures of crystallization, confirming the cold, pristine nature of these lines of sight. After the 
dominant species are removed, there are residuals that suggest the presence of minor species such as 
HCOOH and possibly NH3. Clearly, models of star formation should begin with dust models already 
coated with a fairly complex mixture of ices. 
Subject headings: ISM: molecules, astrochemistry, infrared 



1. INTRODUCTION 

Infrared absorption studies of protostars embedded in 
dense clouds have shown that dust grains along these 
lines of sight have icy mantles (e.g., Willner et al. 1982, 
Tielens et al. 1984, Allamandola et al. 1992, Whittet et 
al. 1996, Boogert et al. 2004). Heating by the protostar 
and energetic photons can affect the ice composition by 
sublimation and by triggering chemical reactions (e.g., 
Gerakines et al. 1996, Ehrenfreund & Charnley 2000, 
Schuttc & Khanna 2003). Knowledge of the ice composi- 
tion in quiescent dense clouds is required in determining 
the amount of processing ices undergo during star forma- 
tion. This 'baseline' can be obtained by observations of 
field stars lying behind molecular clouds. These obser- 
vations also help constrain models of chemical evolution 
during the star forming process (e.g., Lee et al. 2004). 

The 3 /im absorption band of H2O ice was observed 
toward stars behind the Serpens (Eiroa & Hodapp 1989) 
and Taurus dark clouds (Whittet et al. 1988, Smith et al. 
1993, Murakawa et al. 2000). These studies indicate that 
H 2 ice is formed deep in the clouds, at visual extinctions 
Ay > 3 magnitudes (the ice formation threshold). Solid 
CO was observed toward background stars as well (e.g., 
Whittet et al. 1985, Chiar et al. 1994, 1995). Its forma- 
tion threshold is significantly larger (Ay—6-lb mag) due 
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to the lower sublimation temperature of solid CO. 

Studies of ices toward background stars have been lim- 
ited to bands below 5 /im because of telluric absorption 
and the fact that stellar fluxes drop rapidly with in- 
creasing wavelength. Observations of background stars 
have become increasingly feasible with the Infrared Space 
Observatory (ISO) mission and, in particular, with the 
launch of the Spitzer Space Telescope (Werner et al. 
2004). With ISO, solid C0 2 at 4.25 (im was observed to- 
ward two Taurus background stars (Whittet et al. 1998; 
Nummelin et al. 2001), indicating that radiation from 
nearby protostars is not required to form this species. 
Recent observations with Spitzer detected the CO2 bend- 
ing mode at 15 /im toward background stars (Bergin et 
al. 2005). Here we present observations of ices toward 
three background stars over the full 5-20 /im range taken 
with the Infrared Spectrograph (IRS; Houck ct al. 2004) 
aboard Spitzer. We assess the complete ice inventory in 
quiescent clouds and compare it to observations toward 
protostars. 

2. OBSERVATIONS AND REDUCTION 

The background stars discussed here include a source 
behind the Serpens dark cloud, CK 2, and two sources 
behind the Taurus dark cloud, Elias 13 and Elias 16. The 
observations are part of the "c2d" legacy program (Evans 
et al. 2003). All sources were observed with the short 
wavelength, low resolution module (SL; A = 5-14 /xm; 
R = A/AA = 64-128), with on-source integration times 
of 28 sec per spectral order. CK 2 and Elias 13 were 
also observed with the short wavelength, high resolution 
module (SH; A = 10-20 /im; R — 600), with integration 
times of 240 sec and 60 sec per spectral order, respec- 
tively. Spectra of Elias 13, Elias 16 and CK 2 were part of 
AOR# 0005636864, 0005637632, and 0011828224. The 
SH spectrum of Elias 16 has been published by Bergin 
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Fig. 1. — The spectra of Elias 13, Elias 16 and CK 2 with the ex- 
tincted photospheres (dashed lines) used for normalization. Near- 
infrared complementary data are shown in triangles, diamonds and 
plus signs, respectively. The spectra are scaled along the flux axis 
with the numbers given in brackets. 



et al. (2005) and was part of AOR# 0003868160. The 
data were reduced using the Spitzer Science Center (SSC) 
pipeline version SI 1.0. 2 (S12 for SH in Elias 16) to pro- 
duce the 2-dimensional Basic Calibrated Data. Subse- 
quently, customized source extractions were performed, 
including the subtraction of extended background emis- 
sion. The spectra were then defringed with sine-wave 
fitting routines (Lahuis & Boogert 2003). 

Figure ^ shows the Spitzer spectra of the observed 
background stars, complemented by near-infrared (NIR) 
broad band photometry (2MASS 7 ). In addition, for CK 
2, Spitzer IRAC (S. T. Megeath, priv. comm.) and 
ground-based L' band photometry (Churchwell & Koorn- 
neef 1986) are included and for Elias 16 the 2-5 /jm ISO 
spectrum is shown (Whittet et al. 1998). 

In order to put the data on an optical depth scale and 
analyze the ice and dust features, each spectrum is nor- 
malized to the spectrum of an extincted late-type giant 
taken from the ISO database (Sloan et al. 2003). A black- 
body is used at wavelengths below 2.5 /zm to fit the NIR 
photometry. The extinction law used for Serpens is of 
the form A\ ~ A -19 (Kaas et al. 2004), whereas that for 
Taurus has a shallower power law, A\ ~ A -1,7 (Whittet 
et al. 1988). Indebetouw et al. (2005) find a shallower 
slope on the extinction curve beyond 6 ^m. This flat- 
tening is partly due to the silicate and ice features which 
we account for separately and the powerlaw is a good 
approximation to the extinction at longer wavelengths. 
A\ is taken to be 9.1 x (Rieke & Lebofsky 1985). 
Spectral types K4 III, G9 III, and K3 III were adopted 
for CK 2, Elias 13 and Elias 16, respectively. These val- 
ues are within the range of types given in Chiar et al. 
(1994) for CK 2 and close to the K2 III type assigned 
by Smith et al. (1993) for Elias 13 and 16. While the 
selected spectral types give the best fit (in removing the 
photospheric CO and SiO bands at 5 and 8 /im), they are 
uncertain by a few subclasses resulting in an abundance 
error of 20% on the strong features, especially since the 
available ISO database has limited coverage. 

3. RESULTS 

7 This publication makes use of data from the Two Micron All 
Sky Survey, which is a joint project of the University of Mas- 
sachusetts and IPAC/Caltech, funded by NASA and NSF. 
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Fig. 2. — Background star spectra on optical depth scale (black). 
A fit (dashed grey) of small spherical silicate grains is shown to 
which are added laboratory spectra of H2O and CO2 (at 15 fim) 
at intensities corresponding to the column densities mentioned in 
Table 1. 
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Fig. 3.— CO2 feature for a) CK 2, b) Elias 16 and c) Elias 13 
compared to laboratory spectra composed of H2O:CO2=10:l (dot) 
and 1:1 (dash) at 10 K and CO:N 2 :CO 2 =100:50:20 (thin solid) at 
30 K. The grey line shows the best composite fit. Note that the 
vertical offset is due to the H2O libration mode. 



3.1. C0 2 and H 2 Ices 

The CO2 15 /im feature is very strong toward CK 2 
and Elias 16, and weakly detected toward Elias 13. The 
abundance of CO2 relative to H2O is 33% toward CK 2, 
higher than the ^20% abundance seen toward the Taurus 
sources. The derived CO2 column densities toward Elias 
13 and Elias 16 agree within errors with those obtained 
from the 4.25 /im feature (Nummelin et al. 2001; Table 
1). The bottom of the 15 /im CO2 feature appears single 
peaked toward all sources (Fig. [3J and does not show 
double dips due to crystallization as some protostars do 
(Gerakines et al. 1999, Boogert et al. 2004). The pro- 
file of this band toward Elias 16 is fitted in Bergin et al. 
(2005) with the sum of the polar (H 2 0:C0 2 =7:1) and 
the apolar (CO:C02=4:l) mixtures. The polar, H2O- 
rich mixture accounts for 85% of the CO2 column den- 
sity. However, the H2O column density assumed in this 
fit overestimates the observed value by 30%. We require 
that both the observed H2O and CO2 column densities, 
as well as the CO2 band profile are matched. We use 
a combination of two polar mixtures H20:C02=1:1 and 
H2O:CO2=10:l (with the ratios of the two mixtures: 2:1 
for CK 2, 1.3:1 for Elias 16 and 1:0 for Elias 13) and the 
apolar mixture CO:N 2 :CO 2 =100:50:20, all at low tem- 
perature (Ehrcnfreund et al. 1997). Satisfactory fits 
are obtained for polar fractions of 78% in CK 2, 84% 
in Elias 16, and 87% in Elias 13, comparable to the frac- 
tions found by Bergin et al (2005) despite the different 
mixtures used. 

Using the H2O ice column densities from the 3 /im 
band (Table 1) and laboratory spectra of pure H2O ice 
(Hudgins et al. 1993), the H2O bending mode contributes 
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Fig. 4. — a)The spectrum of CK 2 after removing silicate ab- 
sorption is shown in black. The dash grey line is the pure H2O- 
ice contribution to the 6.0 /im feature while the solid grey is the 
H20:C02 mixture, b-d) The dotted lines show the spectra after 
removing the silicate contribution. The solid lines show the spectra 
after removing silicate and H2O contributions, b) For CK 2, the 
dash grey line shows a fit to HCOOH (5.85 and 8.2 /im), NH+ (6.85 
(im) and NH3 (6 /mi), c) For Elias 16, the dash grey line shows a 
fit to NH^~ and NH3. d) For Elias 13, no remaining features are 
found. 



to 77% and 69% of the observed 6.0 /jm absorption fea- 
ture for Elias 16 and CK 2 (Fig. EJ) . However, the peak 
position, width, and strength of this band change signif- 
icantly when H2O is diluted. For example, compared to 
pure H2O, the mixture H20:C02=1:1 shifts the peak to 
longer wavelengths by ~0.1 (tm, and increases the peak 
by a factor of 2.4, but the extensive long wavelength wing 
remains unchanged. Using mixtures that fit the 15 /jm 
C0 2 band (Fig. |3J, ~85-100% of the 6.0 /mi feature 
can be explained by H2O (Fig. 0J. The strong libration 
mode of H2O explains much of the excess absorption in 
the 12-13 /urn region, but due to severe blending with 
the silicate absorption feature, residuals in that spectral 
region are hard to interpret (§3.3). A model of astro- 
nomical silicates (Weingartner & Draine 2001) is used to 
fit the 10 and 20 /jm features for all sources (Fig. [2J, but 
this silicate model may not be unique. 

3.2. The 6.85 /im Band and Other Ices 

A strong feature at 6.85 /tin is detected toward CK 
2 and Elias 16 but not toward Elias 13. This is the 
first detection of this band toward background stars. It 
is commonly observed toward protostars (Keane et al. 
2001) and often attributed to the NH|" ion (see Schutte 
& Khanna 2003; §4). Regardless of the identification, 
the profile may be a powerful tracer of the thermal his- 
tory of the ices. Figure [5] shows the decomposition of 
the 6.85 /jm feature into the short and long wavelength 
components used by Keane et al. (2001) to characterize 
this band. The phenomenological separation is meant 
to represent NH^ bands at different temperatures. The 
ratio of the peak optical depths of the two components 
(short/long) is 2 for Elias 16 and 1.2 for CK 2. If instead 
the equivalent widths are compared, the ratios are 1.6 
and 0.9, respectively. In either case, the profiles resem- 
ble those of the protostars with the coldest sight-lines 
(e.g., NGC 7538 IRS9). Its equivalent width is used to 
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Fig. 5. — The decomposition of the 6.85 /mi feature into short 
(dot) and long wavelength (dash) components (Keane et al. 2001) 
for a) Elias 16 and b) CK 2. The solid grey line indicates the sum 
of the short and long wavelength components. The solid black 
line represents the spectrum after subtraction of silicate and H2O 
contributions. 



TABLE 1 
Column Densities and Abundances 



Species 


Unit 


CK 2 


EL 13 


EL 16 


HH46 1 


B5 1 


CO 


%H a O 


36 2 , 


9 4 


26 4 


20 


43 






57 3 










CO2 


%H 2 


33 


<15, 


18 5 ,22 6 , 


.32 


37 








22 s 


19 7 ,24 






HCOOH 


%H 2 


1.9 






<8.7 


9.3 


CH3OH 


%H 2 


<2.1 


<0.8 


<2.3 


7.0 


< 2.3 


NH 3 


%H a O 


<8 


<6 


<8 


17 


<9.1 


NH+ 


%H 2 


10.8 




5.5 


9.6 


<12.7 


CH 4 


%H 2 


<3 




<3 


4 




OCN- 


%H 2 






<2.3 8 


<0.7 


<0.5 


H 2 


10 18 cm^ 2 


3.5 9 


1.0 10 


2.5 10 


8.4 


2.6 


H 


10 22 cm -2 


6.4 


1.9 


3.6 


5.0 


3.2 



References. — 1. HH46 IRS and B5 IRS1 from Boogert et al. 
2004, and in preparation; 2. Chiar et al. 1994; 3. Pontoppidan, et 
al. 2003; 4. Chiar et al. 1995; 5. Whittet et al. 1998, 6. Nummelin 
et al. 2001; 7. Bergin et al. 2005, 8. Whittet et al. 2001; 9. Eiroa 
& Hodapp 1989; 10. Whittet et al. 1988, Smith et al. 1993 



calculate the column density for NH|~ (Table 1) using the 
band strength from Schutte & Khanna (2003). 

3.3. Other Ices and Hydrogen Column Densities 

H2O ice accounts for up to 85% of the 6.0 /im absorp- 
tion band for CK 2 (Fig. §3.1). Some excess ab- 
sorption remains, in particular on the short wavelength 
side. This is observed towards protostars as well, and 
attributed in part to absorption by HCOOH (Schutte 
et al. 1999, Keane et al. 2001). The 5.85 /mi band is 
the strongest HCOOH band in our observed wavelength 
range, that at 8.2 /tm is slightly weaker. For CK 2, the 
laboratory spectrum of pure solid HCOOH fits the 5.85 
/im feature with a discrepancy at 8.2 /im, which may indi- 
cate over-correction of the photospheric SiO band (§2). 
Table 1 shows upper limits and tentative detections of 
weak features in the 7-13 /im spectral region. 

The uncertainty in the shape of the silicate feature is 
large in some places (e.g., 20% at 10-13 /im) and hard 
to quantify in others. For ice abundance determinations, 
the hydrogen column density is usually calculated from 
the relation A H « 1.87 x 10 21 x Ay cm" 2 for R v = 3.1 
(Draine 2003). As outlined in §2, by fitting ISO tern- 
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plate spectra, extinctions Ay of 34™, 19™, and 10™ are 
derived for CK 2, Elias 16, and Elias 13 respectively. 
The observed peak optical depth of the silicate band and 
the relation Ay /tq,t= 18.5 ±2 (Draine 2003) give similar 
values for Ay. The Ah values listed in Table 1 are calcu- 
lated using the Ay values derived from the photospheric 
fits. 

4. CONCLUSIONS 

Our Spitzer spectra of background stars show clear de- 
tections of absorption features at 6.0 and 6.85 /im that 
had previously only been seen toward YSOs, imposing 
new constraints on the origin of the 6.85 /jm band. The 
strength of the 6.85 /im band, scaled to H2O, is similar 
to that seen toward YSOs, as is the factor of 2 vari- 
ation between sight-lines (Table 1; Schutte & Khanna 

2003) . In one scenario, the 6.85 /im band is explained 
by NH^~ produced by acid-base reactions in ice mixtures 
containing NH3 and HNCO. In laboratory experiments 
such reactions occur at temperatures as low as 10 K with 
conversion factors between 15% and 100% depending on 
ice mixtures and temperature (van Brockhuizen et al. 

2004) . The strength of the 6.85 /im band is a factor of 2 
larger toward CK 2 compared to Taurus. This sight line 
probes a very cold region, as evidenced by the very high 
abundance of the volatile apolar CO ice (T subl ~20 K) 
as well as the smooth profile of the 15 /im CO2 band. In 
contrast, the line of sight of the YSO HH46 IRS shows 
evidence for the ices to have undergone thermal process- 
ing (Boogert et al. 2004), but the 6.85 /im band is not 
unusually deep (Table 1). Thus, if the 6.85 /im band is 
due to NH^, variables other than temperature, such as 
the initial HNCO and NH3 abundances, must play roles 
in determining its strength. We stress that the identifi- 
cation of the 6.85 /im band with NELj" is tentative and 
more evidence, including correlation with the bands of 
counter ions such as OCN~ (van Brockhuizen et al. 2005) 
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or HCOO~ is required. 

Most (^75%) of the 6.0 /im band toward background 
stars is explained by pure H2O ice, comparable to the 
percentage toward most massive YSOs (Keane et al. 
2001). While H2O mixed with CO2 can explain the re- 
maining absorption in the Taurus sources, the residual 
toward CK 2 has a peak absorption wavelength consis- 
tent with HCOOH (Fig. Its abundance would be a 
few % of H2O, comparable to that seen in several high- 
mass YSOs (Keane et al. 2001) but not as high as seen 
in some low- mass YSOs (Boogert et al., in preparation). 

Dust grains have accumulated rather complex icy man- 
tles in opaque regions of molecular clouds before star for- 
mation begins, a point which must be included in models 
of star formation. Also, the effect of freeze-out on the 
thermal balance has been studied by Goldsmith (2001). 
For the two stars with extinctions above 15 mag, the 
abundances relative to H2O ice are within the range seen 
toward embedded objects. From the abundances in Ta- 
ble 1, the percentages of nitrogen, oxygen and carbon 
locked in ices for Elias 16 and CK 2 are 35-37% N, 28- 
30% O and ~12% C. Further work on larger samples of 
background stars will elucidate the dependence of the 
ice composition on cloud conditions and history. Such 
surveys are now possible with Spitzer/IRS for sources as 
weak as 10 mJy (9.5 mag) at 8 /im and Ay of up to 50 
mag. 

The authors thank S. T. Megcath and H. Eraser for 
input. Support for this work, part of the Spitzer Legacy 
Science Program, was provided by NASA through con- 
tracts 1224608, 1230779, 1256316, issued by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under NASA contract 1407. Astrochemistry in 
Leiden is supported by a NWO Spinoza and NOVA 
grant, and by the EU RTN-PLANETS (HPRN-CT-2002- 
00308). 



Nummelin, A., Whittet, D. C. B., Gibb, E. L., Gerakines, P. A., & 

Chiar, J. E. 2001, ApJ, 558, 185 
Pontoppidan, K. M., et al. 2003, A&A, 408, 981 
Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618R 
Schutte, W. A., et al. 1999, A&A, 343, 966 
Schutte, W. A. & Khanna, R. K. 2003, A&A, 398, 1049 
Sloan, G. C, Kraemer, K. E., Price, S. D., & Shipman, R. F. 2003, 

ApJS, 147, 379 

Smith, R. G., Sellgren, K., & Brooke, T. Y. 1993, MNRAS, 263, 
749 

Tielens, A. G. G. M., Allamandola, L. J., Bregman, J., Goebel, J., 
Witteborn, F. C, & Dhendecourt, L. B. 1984, ApJ, 287, 697 

van Broekhuizen, F. A., Keane, J. V., & Schutte, W. A. 2004, A&A, 
414, 425 

van Brockhuizen, F. A., Pontoppidan, K. M., Fraser, H. J., & van 

Dishoeck, E. F. 2005, A&A, 441, 249 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Werner, M. W. et al. 2004, ApJS, 154, 1 

Whittet, D. C. B., McFadzean, A. D., & Longmore, A. J. 1985, 

MNRAS, 216, 45P 
Whittet, D. C. B., Bode, M. F., Longmore, A. J., Adamson, A. J., 

McFadzean, A. D., Aitken, D. K., & Roche, P. F. 1988, MNRAS, 

233, 321 

Whittet, D. C. B., et al. 1996, A&A, 315, L357 
Whittet, D. C. B., et al. 1998, ApJ, 498, L159 

Whittet, D. C. B., Pendleton, Y. J., Gibb, E. L., Boogert, A. C. A., 

Chiar, J. E., & Nummelin, A. 2001, ApJ, 550, 793 
Willner, S. P., et al. 1982, ApJ, 253, 174 



